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Abstract. This paper presents a method of estimation of 
emission altitudes using observational data - precise mea- 
surements of pulse profile widths at low intensity level. 
The analysis of emission altitudes obtained using this 
method for a large number of pulsars gives constraints 
that should be useful for theory of coherent pulsar emis- 
sion. It seems that radio emission originates at altitudes 
of about few percent of the light cylinder and that they 
depend on frequency v, pulsar period P and period deriva- 
tive P. 



1. Introduction 

It is thirty five years since the first pulsar was discovered, 
but the mechanism of coherent radio emission in pulsars is 
still not well known. In the polar cap model (e.g. Goldrc- 
ich & Julian 1969, hereafter GJ; Rudcrman & Sutherland 
1975), the observed emission of radio pulsars is a coher- 
ent radiation of relativistic particles flowing along dipolar 
field lines. This radiation is generated not far from the sur- 
face of the neutron star (NS). The emission altitude r em , 
which is the distance from NS to a place where radiation is 
generated, is a crucial parameter of emission region. It is 
commonly believed that r em dependes on the observed fre- 
quency. This view, supported by the systematic narrowing 
of average pulse width with increasing frequency (Thors- 
set 1991), is described by a radius-to-frequency mapping 
(RFM). It is obvious that all observational limits for emis- 
sion altitude can be crucial for understading the physical 
mechanism generating pulsar radiation. One of the ways 
to estimate the emission altitude (Cordes 1978, 1993) in- 
volves an analysis of observational data concerning the 
average pulse profiles. This method was adopted and de- 
veloped by Kijak & Gil (1997, 1998, hereafter Papers I 
and II) and Kijak (2001, hereafter Paper III), using pre- 
cise measurements of pulse profile width for the largest 
available data set (about 40 PSRs). Three straightfor- 
ward assumptions were used: (1) the pulsar radiation is 
narrow-band, with a radius-to-frcquency mapping operat- 
ing in the emission region (that is, a narrow-band frequen- 
cies Ais <C Ai/ to t are emitted at given altitude r cm , where 



A^tot ~ 30GHz is a total band of pulsar radiation), (2) the 
elementary pulsar emission is relativistically beamed tan- 
gentially to dipolar magnetic field lines (beaming angle 
9 ~ I/7; 7 ~ 100), (3) the extreme profile wings originate 
at or near the last open dipolar field lines. The estima- 
tion of emission altitude is correct (i.e. valid for describ- 
ing the actual physical phenomenon) as long as the above 
assumptions are true. In other words, if we question one 
of these assumptions, our estimation of emission altitude 
can be questioned too. The first two assumptions are in 
agreement with basic pulsar electrodynamics (GJ; Stur- 
rock 1971; Ruderman & Sutherland 1975) and with the 
observed narrowing of profiles with increasing frequency 
(Thorsset 1991; Mitra & Rankin 2002). These two assump- 
tions are commonly accepted. In the third assumption, the 
edge of pulsar radio beam is determined by the boundary 
of the GJ polar cap. This assumption is not widely ac- 
cepted and will be discussed details in §3. 

2. Geometry of emission region 

The pulsar beam originates within the open dipolar field 
lines at altitudes much smaller than the radius of the light- 
cylinder tlc = c/f2 (Lyne & Smith 1998). The last open 
field lines define the polar cap on the neutron star surface. 
The radiating plasma streams outwards along the open 
field lines. Let us consider the case of a neutron star with 
aligned rotation and magnetic axes (see Fig. 1). The equa- 
tion of dipolar magnetic field lines in polar coordinates 
(r, 1?) has a form r = (R 3 /d 2 )-smi) 2 . For the last open field 
line = 7r/2 and r = c/fl, therefore d = ^/R 3 (2ir/cP) and 
the radius of polar cap is r p = 1.45 ■ 10 4 P _1//2 [cm], for the 
neutron star radius R = 10 6 cm. 

2.1. Geometry on the celestial sphere 

The geometry of pulsar radio emission as defined by 
Manchester & Taylor (1977) is presented in Fig. 2. The 
rotation and magnetic axes are inclined by an angle a. 
The observer's angle £ = a + f3, where j3 is the impact 
angle of the closest approach of the line-of-sight (1-of-s) 
to the magnetic axis. The beam radius (corresponding to 
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Fig. 1. Geometry of dipole field lines of the neutron star 



a pulse profile) is p > (3 and the observed pulse width is 
W = 2if, where ip is the longitude of the beam edge. The 
fiducial plane containing the rotation f2 and magnetic m 
axes defines the so called fiducial phase, which is usually 
taken as the longitude <p — 0° (Fig. 2). 



2.2. The method of estimating emission altitude 

The opening angle of the pulsar beam corresponding to 
the pulse width W is given by 

1/2 



p = 2 sin 



sin 



sin a • sin(a + j3) + sin 2 — 



W 
4 



(1) 



where a is the inclination angle between the rotation and 
magnetic axes and /3 is the impact angle (Gil, Gronkowski 
& Rudnicki 1984). This formula is independent of the 
shape of the pulsar beam (circular, elliptical, patchy, etc.) 
when the beam is symmetric with respect to the fiducial 
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Fig. 2. Geometry on the celestial sphere. 

plane. This of course assumes that the edge of beam is 
close to circular or elliptical in shape. 

The opening angle p expressed above is also the angle 
between m axis and the tangent to magnetic field lines 
at points where the emission corresponding to the appar- 
ent pulse width W originates. For dipolar field lines (see 
Fig. 1) 

p= l.2A° s(r/Rf 2 P- 1 ' 2 , (2) 

where P is the pulsar period, r is the emission altitude 
(at which radiation at a given frequency v is generated), 
R = 10 6 cm is the neutron star radius and the mapping 
parameter s = d/r p describes the locus of corresponding 
field lines on the polar cap (s = at the pole and s = 1 
at the edge of the polar cap. Finally, we calculate the 
emission altitude as follows 




3. Precise measurements of pulse profile and the 
edge of the radiation beam 

The most important question is what we mean by the 
edge of pulsar beam and the edge of an average pulse pro- 
file. First of all, the edge of pulsar beam is related to the 
edge of a pulse profile in the following term: when the pul- 
sar beam starts passing through the radiotelescope, then 
we start observing a weak coherent signal in the receiver. 
The moment at which the pulsar beam enters the radiote- 
lescope corresponds to the moment at which the recorded 
noise rapidly becames a coherent signal (and signal be- 
camcs a noise when the beam leaves the radiotelescope). 



In studies on emission regions (Papers I, II and III) we 
used precise measurements of the longitude correspond- 
ing to these events. The egde of the beam is defined as 
the longitude at which recorded signals are distinguished 
from the noise and have the instensity at least 0.1% of the 
maximum intensity. Since the pulsar radiation is relativis- 
tically beamed along dipolar field lines, we can attempt 
to calculate which field lines should be tagged as coming 
from the emission region. This should be easiest for the 
last open magnetic field lines, which arc believed to be as- 
sociated with the lowest detectable level of radio emission 
- i.e. at the profile wings. 

The radio emission from pulsars is characterized by 
a short time of duration, typically 3-10% of its period. 
However, a few objects demonstrate the emission through 
the entire pulsar period. Hankins et al. (1993) showed, 
using the VLA, that the emission through the entire pulsar 
period is observed only in cases when the inclination angle 
a is smaller than the beam radius p (see Fig. 2), that is, 
the 1-of-s is always inside the emission beam (Gil 1985). 
In general, there is a transition from the noise to the weak 
emission at the extreme profile wings (Fig. 3 and 4). In 
our analysis of emission altitudes (Paper I, II and III), 
we measured pulse widths at the extremely low intensity 
level corresponding to about 0.05 per cent of the maximum 
intensity (sec Fig. 3), using the polar- log-scale technique 
(Hankins & Fowler 1986). 

It is not certain that the entire GJ polar cap is involved 
in the emission process thus, the edge of pulsar beam is 
emitted at or near the last open field lines (which corre- 
spond to s = 1). It is possible that a forbidden area exists 
at the outer parts of the GJ polar cap. Evidently, we can- 
not give a theoretical answer to these questions. Recently, 
Gangadhara & Gupta (2001) found new outermost weak 
components in the pulse window of pulsar B0329+54 at 
frequencies of 320 and 610 MHz. The separation A(j> be- 
tween the newly detected components is 34° . In Fig. 5 the 
profile of PSR B0329+54 at 320 MHz is presented, which 
was used by Kijak & Gil (1998) in their measurement of 
pulse width and corresponding estimation of emission al- 
titude (Paper II, Fig. 2). Let us note that the pulse width 
measurement W at 0.05% level is about 49°. A compari- 
son of these two measurements (A<j> = 34° and W = 49°) 
clearly shows that the profile width presented in Fig. 5 
includes new components found by Gangadhara & Gupta 
(2001). Thus, the reported new outermost components in 
the profile of PSR B0329+54 do not contradict our as- 
sumption that s « 1 at the pulse edges. Our method of 
profile width measurements has proven to be sensitive to 
very weak radio emission at profile wings originating at 
the edge of the beam, most probably close to the last 
open field lines. These outermost components are clearly 
visible in Fig. 5 at the level of about 0.1% of the maximum 
intensity. 
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Fig. 3. (a) A plot of profile of PSR B0329+54 at 1.4 GHz, 
(b) the profile of the same pulsar in the log-scale and the 
polar coordinates. The dynamic range is 10 4 . 

4. Emission altitude estimation 

In our research of radio emission regions we used the high 
quality average pulse profiles in the wide frequency range. 
We selected pulsars for which some information about a 
and (3 angles was available. For this purpose, observations 
were performed between 1995 and 1998 using the Effels- 
berg 100-m radio telescope (Lorimer et al. 1998) and the 
Jodrell Bank data (Gould & Lyne 1998) were used as well. 
As a result, the high quality average profiles were obtained 
with the signal-to- noise ratio S/N~ 1000 in the frequency 
range between 0.3 and 20 GHz. These profiles were used 
to obtain precise measurements of pulse profile widths and 
to estimate the emission altitudes. 

4-1- Results 

The result of emission altitude analysis (Papers I, II and 
III) indicated clearly that the emission altitude r cm de- 
pends not only on the observing frequency v, but also 
on the basic pulsar parameters: pulsar period P and pe- 
riod derivative P. Here, we present the analysis of cmis- 



Fig. 4. The syntetic data with pulsar signal which comes 
from open field lines. The noise was obtained using a ran- 
dom generator. Error box is marked as an example (see 
details in Paper I). 



sion altitudes at 1.4 GHz for the largest available data 
set (37 PSRs). An apparent period dependence of emis- 
sion altitudes is visible in Fig. 6, where a formal fit to all 
data points gives r cm = (45 ± 3)i?P°- 37±0 05 . This result 
confirms the previous analysis of the emission altitude at 
different frequencies for a smaller data set presented in 
Papers I and II. Figure 6 also shows a weak dependence 
r om on pulsar characteristic age tq (in 10 6 yrs), which is 
discussed in Paper III (sec its Fig. 2 and Table 2). For this 
analysis, a formal fit to the same data set (37 PSRs) at 



1.4 GHz yields r cm = (33 ± l)i?r 6 



-0.07±0.03 



. In Paper II, 



the RFM was derived for 16 pulsars covering the frequency 
range between 0.3 and 20 GHz. The average frequency de- 
pendence on emission altitude is represented by the formal 
fit r em cx jy- - 26 ±o.09^ Generally, the pulsar emission alti- 
tude depends on frequency v 1 pulsar priod P and pulsar 
age r. 
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Fig. 6. Emission altitude at 1.4 GHz versus pulsar period. 
Open circles correspond to young objects (see Paper III 
for details). 



Fig. 5. A pulsar profile of B0329+54 at 320 MHz in the 
log-scale and the polar coordinates. The data come from 
Jodrell Bank Observatory. 

4-2. The semi- empirical formula for emission altitude 

The radial location of the emission regions was discussed 
in Papers I, II and III. A general form of the formula for 
the emission altitude is 

r KG =AR^ z r 6 b P c , (4) 

where A *~ 55, R = 10 6 cm, vqyiz is observing frequency 
(in GHz), T6 is characteristic age (in 10 6 yrs), P is pulsar 
period (in sec). 

The study on emission region in Paper I carried out for 
about 10 PSRs at two frequencies v = 0.4 and 1.4 GHz 
yielded A = 55 ± 5, a = 0.21 ± 0.07, b = 0.07 ± 0.03 and 
c = 0.33 ± 0.05. Here we verify these values for a larger 
data set, as well as for a wider frequency range. Taking 
into account results of § 4.1, we obtain 

_ frn 4- 1 n\ p„-0.26±0.09 -0.07±0.03 p0.37±0.05 /r\ 

r-KG - (50 ± W)Kv Gn!e t 6 P , (5) 

or 

r KG « (40 ± 8 ) i?I/ -0.26±0.09p0.07±0.03p0.30±0.05 ; (6) 

where rg = 16P/P_i5 an d P-w = -P/10 15 is the period 
derivative. The uncertainties given for various quantities 
in the above equations follow from systematic errors in 
the pulse width W measurements and random errors in 
estimations of the inclination angle a and the impact an- 
gle [3. 



5. Discussion and Conclusions 

As it follows from Eqs. (5) or (6), the radio emission of pul- 
sars originates at the narrow altitude range below 10 per 
cent of the light-cylinder radius. In millisecond pulsars, 
the radio emission region is much more compact, consis- 
tent with their small light-cylinder radii, and is located 
closer to the neutron star surface (see Paper II). Previ- 
ously, it was thought that millisecond pulsars have differ- 
ent characteristics of radio emission from those of normal 
pulsars. However, it was recently argued (Gil & Sendyk 
2000; Kramer 1999, and this proceedings) that millisec- 
ond pulsar emission properties do not differ from those of 
typical pulsars. 

Recently, Gil et al. (2002) analysed simultaneous dual- 
frequency single pulse observations of PSR B0329+54 at 
240 and 610 MHz. The phase shifts of the leading and the 
trailing conal components are not equal at these frequen- 
cies. This is caused by the retardation-aberration effects. 
In fact, since the conal beams at different frcquenics arc 
emitted at different altitudes, lower frequencies arrive ear- 
lier. The retardation-aberration shift of 0.39° translates 
into a difference of emission altitudes of Ar ~ 2.3 x 10 7 cm. 
This is consistent with the result obtained using the semi- 
empirical formula Arxc = 2x 10 7 cm (Eqs. 5 or 6). 

Mitra & Rankin (2002) found the RFM using outer- 
conal components data. In Fig. 7 we present a comparison 
between their results and calculation from Eqs. (5) or (6). 
Emission altitudes represented by open and filled circles 
were calculated from r cm = h ■ s -2 , where h was the emis- 
sion height for conal component and s was taken as the 
average value from Mitra & Rankin's paper for each pul- 
sar, respectively. The open circles represent calculations 
for the first group of their fitted procedure, and filled cir- 
cles were calculated using the index of —2/3 (see Mitra 
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Fig. 7. The multifrequency study of emission altitude. Data points represent calculations from the paper of Mitra & 
Rankin (2002). The dashed lines show upper and lower limits, calculated from equation (5), corresponding to limiting 
values of parameters in the semi-empirical relationship. 



& Rankin 2002, their Table 4). Figure 7 shows that the 
semi-empirical formula describes quite fr 

well the radio emission region in pulsars, consistenly 
with other results. 

The main conclusions from the work on radio emission 
regions are as follows: 

(i) A radius-to-frcquency mapping operates in the emis- 
sion region. 

(ii) Pulsar radio emission is typically generated at alti- 
tudes lower than a few per cent of the light-cylinder 
radius tlc and the ratio fKG/fLC decreases with in- 
crcsing period. The total size of the emission region is 
below 500i? in longer period pulsars and correspond- 
ingly smaller in shorter period pulsars. 

(iii) The emision region in old pulsars is located at corre- 
spondingly lower altitudes than in young pulsars with 
approximately the same period. 

Qq cx 
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